Phonon effects on the radiative recombination of excitons in double quantum dots 
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We study theoretically the radiative recombination of excitons in double quantum dots in the 
presence of carrier-phonon coupling. We show that the phonon-induced pure dephasing effects 
and transitions between the exciton states strongly modify the spontaneous emission process and 
make it sensitive to temperature, which may lead to non-monotonic temperature dependence of the 
time-resolved luminescence. We show also that under specific resonance conditions the biexcitonic 
interband polarization can be coherently transferred to the excitonic one, leading to an extended 
life time of the total coherent polarization, which is reflected in the nonlinear optical spectrum of 
the system. We study the stability of this effect against phonon-induced decoherence. 
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I. INTRODUCTION 

Systems composed of coupled nanostructures attract 
much attention as they not only open the path to new 
applications but also provide new ways of controlling the 
states and dynamics of confined carriers. For instance, 
systems build of two closely located quantum dots pro- 
vide the basis for an implementation of excitonic qubits 
with extended life times 1 -, allow one to effect conditional 
quantum control of carrier states 2 -, and have been pro- 
posed as a source of entangled photons^. Carrier transfer 
in such double quantum dot (DQD) structures can also 
be used for initializing the spin state of a dopant Mn 
atom in one of the dots 4 . 

The optical properties of the DQD structures are much 
more complex than those of a single quantum dot (QD), 
which is manifested in the occupation and coherence dy- 
namics observed in optical experiments^—. This is partly 
due to the richer structure of carrier states^— but also to 
the interplay of a few other factors that govern the evo- 
lution of carriers in these systems after an optical excita- 
tion: collective spontaneous emission (radiative recombi- 
nation of excitons} 1 - 2 -, coupling between the dots 1 ^, as well 
as phonon-induced transitions 1 ^— and dcphasingiSr— . 
In our previous work, we have investigated the effect 
of the first two of these factors 1 - 2 - and their role in the 
linear— and nonlinear— response of DQDs. We have 
shown, in particular, that collective effects in the sponta- 
neous emission of uncoupled DQDs disappear as soon as 
the transition energy mismatch between the dots forming 
the system exceeds the spectral line width of the funda- 
mental transition (on the order of micro-electron- Volts), 
which is well beyond current manufacturing possibilities. 
However, coupling between the dots leads to pronounced 
collective effects even for technologically realistic values 
of the energy mismatch on the order of milli-electron- 
Volts. We have also pointed out that phonon-related 
effects can affect the collective spontaneous emission in 
these system o 23 ! 24 . 

In this paper, we study the effects of phonon-induced 
transitions and dephasing on the kinetics of collective ra- 
diative recombination in DQDs. The interplay of these 



two factors has recently been shown to determine the 
recombination properties of a lateral DQD system with 
direct to indirect exciton tunneling 25 . Here, we con- 
sider a vertically stacked system with realistic values of 
the parameters and focus on the experimentally observ- 
able features in the system dynamics assuming the spa- 
tially direct nature of the excitons (no electric field ap- 
plied). We study the decay of the exciton occupation 
in a DQD in the presence of phonons and show that 
it becomes strongly dependent on temperature due to 
phonon-induced thermalization of single-exciton occupa- 
tions and pure dephasing processes. For negative cou- 
pling (characteristic of tunneling), this leads to an in- 
crease of the exciton lifetime with temperature and to its 
non-monotonic behavior in a typical situation when ther- 
mally activated non-radiative decay processes are present 
at higher temperatures. We show also that for specific 
values of the system parameters, coherent transfer of in- 
terband polarization between the biexciton and exciton 
transitions can take place, which is reflected in the non- 
linear optical spectroscopy of these systems, and we study 
the stability of this process against phonon-induced de- 
coherence. 

The paper is organized as follows. In Sec. [ill we define 
the model under study. Next, in Sec. Mil we define the 
method of numerical simulations and present a prelimi- 
nary discussion of the most important phenomena that 
will determine the evolution, based on analytical solu- 
tions in the limiting cases. Sec. IIVI contains the results 
of our modeling and their discussion. Finally, Sec. fVl 
concludes the paper. 



II. MODEL 

We consider a self- assembled structure composed of 
two vertically stacked QDs with relatively similar sizes in 
the absence of an external electric field. The lowest sector 
of exciton eigenstates is then formed by spatially direct 
states, that is, configurations in which the electron-hole 
pairs reside in the same do t 11 ' 26 . We restrict the discus- 
sion to the lowest exciton state in each dot and assume 
that the spins of the carriers are fixed, corresponding to 
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bright exciton configurations. Under these assumptions, 
each dot can be either empty or occupied by one exci- 
ton and the model includes four basis exciton states: the 
state without excitons, |0), the states |f) and |2) with 
an exciton in the first and second dot, respectively, and 
the "molecular biexciton" state \B) with excitons in both 
dots. The carriers confined in the QDs interact with the 
longitudinal acoustic phonon modes via the deformation 
potential coupling and with the photon (radiative) envi- 
ronment via the dipole coupling. 

The total Hamiltonian of the system is 



H = Hi 



DQD 



H, 



rad 



c— ph 



rad ■ 



The first term describes exciton states in the DQD struc- 
ture and has the form 



(S) 
DQD 



J2 ei\i){i\ + (ei + e2+E B )\B){B\ 



i=l,2 



+V(|1)<2|+|2}<1|), 



(1) 



where t\ and £2 are the fundamental excitonic transi- 
tion energies in the two dots, V is the coupling between 
the dots, which may originate either from the Coulomb 
(Forster) interaction or from tunnel coupling (it can be 
assumed real), Eb is the biexciton shift, and the super- 
script "S" denotes the Schrodinger picture. Electron and 
hole wave functions in the dots 1 and 2 are modeled by 
identical anisotropic Gaussians with identical extensions 
I in the xy plane and l z along the growt axis z for both 
particles, 



ip lj2 (r) ~ exp 



y 



{z±D/2f 



21 2 



21 2 



(2) 



where the upper (lower) sign refers to ip\ (1P2) and D is 
the distance between the dots. 

The phonon modes are described by the free phonon 
Hamiltonian 

fc 

where bk,b\, are the bosonic operators of the phonon 
modes and ujk are the corresponding frequencies. We 
assume a linear dispersion relation for phonons. Inter- 
action of carriers confined in the DQD with phonons is 
modeled by the independent boson Hamiltonian 



Hr.-ph — 



i=l,2 



+ 



\B)(B\)J2fL l \bl 



(3) 



(12) 

where f k ' are system-reservoir coupling constants. For 
Gaussian wave functions as in Eq. ^ty, the coupling con- 
stants for the deformation potential coupling between 
confined charges and longitudinal phonon modes have the 



form f k = /fce ±lfczI) / 2 , where 



/fc = (Cc — C7h)l 



2gvc 



exp 



I ~k y + I k 1 



Here v is the normalization volume, k±/ z are momentum 
components in the xy plane and along the z axis, cr c ^ h 
are deformation potential constants for electrons/holes, 
c is the speed of longitudinal sound, and g is the crys- 
tal density. We assume that off-diagonal carrier-phonon 
couplings are negligible due to small overlap of the wave 
functions confined in different dots. 

The third component in our modeling is the radiative 
reservoir (modes of the electromagnetic field), described 
by the Hamiltonian 



fc.A 



where a*, a, aj, A are photon creation and annihilation op- 
erators and Wk are the corresponding frequencies (A de- 
notes polarizations). The QDs are separated by a dis- 
tance much smaller than the resonant wavelength so that 
the spatial dependence of the EM field may be neglected 
(the Dicke limit). The Hamiltonian describing the inter- 
action of carriers with the EM modes in the dipole and 
rotating wave approximations is 



-rad 



= E_ 



fc,A 



(4) 



with 



and 



12X51 + 10X21 + 11X51 



9k,\ = id ■ e x (k) 



2eo£rV ' 



where d is the interband dipole moment, Eq is the vacuum 
permittivity, e x is the dielectric constant of the semicon- 
ductor and e\(k) is the unit polarization vector of the 
photon mode with the wave vector k and polarization 
A. For wide-gap semiconductors with the band gap on 
the order of an electron- Volt, zero-temperature approx- 
imation may be used for the radiation reservoir at any 
reasonable temperature. 

In our numerical simulations (to be presented in 
Sec. IIV[) . we take the parameters corresponding to a 
self- assembled InAs/GaAs system: er — t7h = 9 eV, 
p = 5350 kg/m 3 , c — 5150 m/s, the wave function param- 
eters I = 4.5 nm, l z — 1 nm, D — 6 nm, and the radiative 
recombination time (for a single dot) 1/Tq = 800 ps. 



III. EVOLUTION 

The evolution of the DQD density matrix in the inter- 
action picture is given by the equation 



P = £rad[p] + £ P h[p], 



(5) 



where p is the density matrix of the exciton subsystem 
and the two terms on the right-hand side account for the 
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radiative and phonon-induced dissipation and dephasing 
effects. Eq. ([5]) is vaiid in the leading order in both the 
couplings, where the two dissipation channels are addi- 
tive. 

The interaction of carries with the radiation reservoir 
leads to a Markovian evolution which can be described 
by a Lindblad generator. For a double dot, the analysis 
based on the Weisskopf-Wigner procedure 1 - shows that 
the correct general description of the collective sponta- 
neous emission is provided by the equation 



where 



P 



(R) 



"DQD' P 



4> (R) ] 



where 



Here the superscript "R" denotes the "rotating frame" 
picture defined by the unitary transformation e lS ^\ with 



S{t) 

where E = (ei 



E 



t(|l)<l| + |2X2|+2|B)< J B|), 



e 2 )/2. That is, 



P 



(R) 



e iS(t) JS) £ -iS(t) 



where p*- s -* denotes the density matrix in the Schrodingcr 
picture, and 



rr(R) _ JS(t) tt(S) iS(t)_ n d 
"DQD — c n DQD e 110 

= A(|l)(l|-|2)(2|) + y(|l)(2| 
+Eb\B)(B\, 



|2)(1|) 



(0) 



where A = (ei — £2)/2 is the parameter describing the 
transition energy mismatch between the dots. Since 
i?DQD] = the density matrix in the interaction 
picture with respect to £?dqd is given by 



p = e iH^ u t/H p (K) e -iH£> D t/h 

and the generator of the dissipative evolution becomes 



S_(t)pS + (t)-i{S + (t)S_(t),p} + 



(7) 



e iH DQn t / h Y,_ e - 



iH^t/h and 



where S_(t) = Y,\_(i) 
T = E 3 \d\ 2 y/e^/(3ireoc 3 h 4 ) is the spontaneous decay 
rate for a single dot. 

For the carrier-phonon interaction, non-Markovian ef- 
fects can be important in some cases. Moreover, there 
seems to be no universal way to extract the Markov 
limit in various physical situations. Therefore, for the 
description of phonon-related effects, we use the time- 
convolutionless equation in the lowest order— 



£ P h[p] = ~ drTr ph [H c ^ ph (t), [i? c _ ph (r) , p(t) 
Jo 



® P P h] 



(8) 



H c - P h(t) 



is the carrier-phonon interaction Hamiltonian in the in- 
teraction picture, p p h is the phonon density matrix at 
the thermal equilibrium, and Tr p h denotes partial trace 
with respect to phonon degrees of freedom. Note that the 
form of -ff c _ph is the same in the interaction and rotating 
frame pictures as [-ff c _ p h, S{t)\ = 0. 

The states |0) and \B) are eigenstates of the DQD 
Hamiltonian, therefore nontrivial evolution generated by 
-^dqd takes place only in the subspace spanned by 
|1), 1 2). From the point of view of the theoretical anal- 
ysis, it is convenient to introduce the parametrization 
A = £ cos 8, V = £sin#, with £ > and 9 € [-7r,7r), 
and to treat £ and as the independent parameters of 
the model. Then, the single-exciton eigenstates of the 
excitonic Hamiltonian can be written as 

|+> = cos-|l)+sin-|2), 

I-) = -sin^|l)+cos^|2), 

and correspond to the energies ±£ . Thus, 2£ is the en- 
ergy splitting between the single-exciton eigenstates of 
H^q D and 6 is the mixing angle of the single-exciton 
states. 

The evolution equation ([5]) can be solved numerically 
at moderate computational cost^. However, some phys- 
ical insight is gained if the two evolution generators given 
by Eq. ([7]) and Eq. ([8]) are treated to some extent ana- 
lytically. 

For the spontaneous emission, we expand the generator 
((Tj) in terms of the eigenstates of H^q B . In the present 
discussion we assume that both the energy splitting pa- 
rameter £ and the biexciton shift Eb are much larger 
than the inverse recombination time. This is reasonable 
as the energy mismatch and the coupling are typically on 
the order of milli-electron- Volts, while the recombination 
time is tr ~ 1 ns, yielding h/r-R, ~ 1 peV. Then, we can 
neglect the oscillating terms proportional to e± l£t / h and 
e ±iE B t/h_ This l eac is to the emission-induced evolution 
of the exciton and biexciton occupations given (in the 
interaction picture) by 

P±± = T± {pbb - P±±) , Pbb = -2T pBB, (9a) 

where p a/ 3 — (a\p\(3) and r± = r (l ± s'mO). The 
optically active coherences (optical polarizations) corre- 
sponding to the exciton and biexciton transitions evolve 
according to 

Po± = ~T ±Po± + T ±P±B e^ 2£ ^ h , (9b) 



p±b = - (r + r±/2) p± B , 



(9c) 



while the evolution of the "dark" exciton-exciton and 
vacuum-biexciton coherences is described by 



Pob 



-FqPob, P±t = ~^oP±t- 



(9d) 



4 



One of the two essential effects following from 
Eqs. (|5a|) - ((9d)) is the appearance of the two different de- 
cay rates T±. As a consequence, one of the single-exciton 
states is brighter (decays faster) and the other one is 
darker (decays slower) than a single exciton state in an in- 
dividual dot (from now on, the terms "dark" and "bright" 
will always refer to the spatial state in the DQD struc- 
ture; the spin configuration is always assumed to corre- 
spond to a bright exciton). For 9 — ±7r/2, one recov- 
ers the perfect superradiance (doubled decay rate) and 
subradiance (no radiative decay) of the two states. The 
former is a precursor of the true superradiance of many 
identical atomic systems^. Note that this limit can be 
reached both for identical dots A — » and for strongly 
coupled ones \V\ 3> |A|, which makes it important also 
in the case of typically strongly inhomogeneous QD sys- 
tems. 

The second effect is related to the oscillating term in 
Eq. (|9bj) . In a general case, it can be expected to have lit- 
tle impact on the secular evolution since the two frequen- 
cies 28 /h and Es/h are typically different. However, the 
biexciton shift depends on the static dipole moment in 
the dots, hence on the spatial arrangement of the carriers, 
and can be tuned e.g. by applying an external axial elec- 
tric field. Therefore, these two frequencies can be tuned 
to resonance and this term becomes non-negligible. Phys- 
ically, this means that the transition from the biexciton 
state to one of the single exciton state has the same en- 
ergy as the transition from the same single exciton state 
to the vacuum state. As follows from Eq. (I9bj) . in such 
case the exciton coherence po± is driven by the biexci- 
ton coherence p±B, hence coherence can be transferred 
from the biexciton transition to the exciton transition. In 
Sec. IIVB| we will show that this transfer leads to a con- 
siderable modification of the decay of optical coherences 
in the signal and affects the nonlinear optical spectrum 
of the system. 

For the phonon contribution, we rewrite Eq. (J3j) in the 
eigenbasis of the exciton states 



a=± k 

-fisin^Ffed+X-l + l-X+l), (10) 

k 

where 



A+) 2 U A 1 ) , • 2 A2) 

fk = cos 2# +Sln 2 f k ' 



/i-^sin^/W+cos^/f, ,11, 



and 



F k = fi 2) -fi 1) =-2if k sm- 



fc.-D 



The coexistence of carrier-phonon couplings diagonal and 
off-diagonal in the carrier eigenstates in this interaction 
Hamiltonian makes the evolution nontrivial 17 . However, 
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FIG. 1: (Color online) The phonon spectral density for a few 
values of the inter-dot distance D. 



the diagonal coupling terms in Eq. (flU)) lead to non- 
Markovian pure dephasing on a picosecond time scale, 
while the off-diagonal terms induce transitions, which 
may be described to a good approximation in the Markov 
limit and are characterized by time scales of several pi- 
coseconds. Using this time scale separation, one can ap- 
proximately treat the two effects separately for the sake 
of interpreting the simulation results to be presented in 
Sec. ED 

The pure dephasing dynamics is due to the lattice re- 
sponse to the ultrafast change of the carrier stat o 29 i 30 . As 
the phonon dynamics depends on the carrier state the lat- 
tice gets correlated to the confined carriers and effectively 
extracts information on the carrier stated. This induces 
dephasing of quantum superpositions of various exciton 
states, which has been studied from the point of view of 
the optical response^ and of the decay of entanglement 
between the dots^. For the present discussion, the es- 
sential effect is the dephasing of superpositions between 
the single exciton eigenstates. Following the approach of 
Refs. Il8ll32i . one finds the evolution of the corresponding 
coherence (neglecting the effect of the off-diagonal terms 
inEq. (HJJJ) 



(+W)H = <+|p(0)|-)e-M*), 



(12) 



where 



h{t) 



4 cos 2 9 



.fk 
Uk 



KD 



h 4 cos 2 



sin 

fk 
Uk 



(1 - cos uj k t)(2n(ijJk) + 1) 



sin 



2 "-IB 



-(2n(u)k) + 1) = h c 



where n(u) are the phonon occupation numbers and t p h 
is the characteristic time of the phonon-induced initial 
dephasing which is of the order of a few picoseconds 18 . 
The factor cos 2 9 reflects the fact that for maximally de- 
localized eigenstates (9 — tt/2) the charge distribution 
is the same for both of them and they cannot be distin- 
guished by the phonon response. 

The real transitions governed by the off-diagonal cou- 
pling terms in Eq. (|10|) in many cases take place on a 
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longer time scale (several picoseconds to tens or hun- 
dreds of picoseconds, depending on the amplitude of 
the coupling) and can be to a good approximation 
described in terms of the Markovian transition rates 
7i_> 3 - = 27rsin 2 OJ(\u)ij\)\n(uij) + 1|, where i,j = ±, 
LUij = (Ej — Ei)/h and the spectral density is defined 
as 

J {(J) = Ifkf sin 2 -^-S(u ~ Wfc). 
fe 

The spectral density for a few values of the inter-dot spac- 
ing D is plotted in Fig. [TJ Apart from the usual cutoff 
at frequencies larger than I/c due to the restriction on 
the momentum non-conservation A(Hk) < h/l, where I 
is the average dot size, it is characterized by the oscilla- 
tions with the period Aw = 2-kc/D which result from the 
double dot structure of the system. These oscillations 
are most pronounced at higher frequencies, where the 
emission along the z axis (the strongest confinement di- 
rection) dominates, which makes the factor sm 2 {k z D /2) 
most efficient. 

The above remarks will provide us with useful guide- 
lines for the interpretation of the simulation results pre- 
sented in the next section. 



IV. RESULTS 

In this section we present and discuss the results of 
numerical simulations of the spontaneous emission form 
double dot system under the influence of phonons. First, 
in Sec. IIV Al we study the decay of DQD occupations. 
Then, in Sec. IIVBI we describe the effect of coherent 
transfer of interband dipole moment and study its sta- 
bility against phonon-induced perturbations. In all the 
simulations, we assume 9 < 0, hence V < and the lower 
single exciton state |— ) is brighter than the higher state 

1+) (T_ > r+). 

A. Occupation decay 

As discussed in Sec. lIIIl the coupling between the dots 
leads to the formation of two delocalized eigenstates with 
different decay rates. For V < 0, the state |+) is darker 
(decays slower) and the state |— ) is brighter (decays 
faster) . This can be seen in Fig. Ufa) (green dashed and 
blue dotted lines), where we show the evolution of the 
exciton occupation in the DQD without carrier-phonon 
coupling (in this case, the evolution depends only on the 
mixing angle 9 but not on £ or the system geometry). 
For a general initial state, both components are present 
in a superposition and the decay is not simply exponen- 
tial, with an average life time in between the two limiting 
cases. An example of such a situation is shown by the red 
solid line in Fig. [2ja) for the particularly interesting ini- 
tial state (|1) + \ 2))/^/2, which corresponds to an optical 
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FIG. 2: (Color online) The decay of the exciton occupation for 
three different initial states as shown for D = 6 nm, 9 = — 7r/6, 
T = 4 K. (a) without carrier-phonon interaction; (b) with 
phonons, for 2£ = 1.77 meV; (c) with phonons, for 2£ = 
3.55 meV. 
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FIG. 3: (Color online) The decay of the exciton occupation 
for the initial state (|1) + |2})/v / 2 for D = 6 nm. (a) dom- 
inated by relaxation (9 = — n/3) with strong carrier-phonon 
coupling (2£ = 1.77 meV); (b) as previously but with much 
weaker phonon effects (2£ = 7.1 meV); (c) dominated by pure 
dephasing (6 = 0,28 = 0). 



excitation by a pulse that resolves the two dots neither 
spatially nor spectrally. 

This kinetics of recombination changes, however, when 
the carrier-phonon interaction is taken into account, as 
shown in Fig. [2]Jb,c). Phonon-induced relaxation redis- 
tributes the occupations on a time scale much shorter 
than the exciton life time, which suppresses the depen- 
dence on the initial state. For strong phonon coupling as 
in Fig. H{b) (corresponding to the first maximum of the 
spectral density, see Fig. [1]) the resulting radiative decay 
is almost insensitive to the initial state. The situation is 
still very similar in the case of energy splitting 2£ cor- 
responding to the first minimum of the spectral density, 
which yields a much weaker effective coupling, as shown 
in Fig. El^c). This shows that the phonon-induced redis- 
tribution of occupations dominates the radiative decay 
over a rather wide range of energy parameters. How- 
ever, for the values of 2£ around the further minima as 
well as beyond the effective cut-off at about 10 meV the 
phonon-related processes become inefficient and the de- 
cay becomes similar to that observed in the absence of 
carrier-phonon coupling. 

The important role of phonons in the radiative recom- 
bination leads to a strong temperature dependence of the 
emission process. This effect is shown in Fig. [3J where 
the exciton occupation decay for the initial optically gen- 
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erated state (|1) + \2))/^/2 at a few temperatures is plot- 
ted for various system parameters. In Fig. [3]^ a) we choose 
9 = — 7r/3 and the value of the energy splitting 2£ close to 
the first maximum of the spectral density [Fig. [Tj , which 
leads to efficient phonon-induced relaxation. The effect 
of the resulting redistribution of occupations is clearly 
seen in the dynamics of recombination: with increasing 
temperature, the occupation of the higher-energy dark 
state increases at the thermal quasi-equilibrium and the 
resulting decay becomes slower. Such a change of the re- 
combination rate with temperature can be observed also 
for higher value of 2£, corresponding to the deep second 
minimum of the spectral density. However, in this case 
the effect is weaker and appears only at higher tempera- 
tures, as a considerable occupation of the higher state ob- 
viously appears only at T ~ £ /fee ■ In Fig. [3Jc) we show 
that a similar effect of extended life time can appear also 
due to pure dephasing between the states localized in the 
two dots^. Here, we set 9 — 0, which excludes occupa- 
tion redistribution. However, the occupation-conserving 
pure dephasing effect remains effective and leads to de- 
struction of the coherence between the two single exciton 
states according to Eq. (TL2")) , where, in the present case, 
|+) = |1) and |— ) = 1 2). Due to this dephasing, the 
initially bright state is partly turned into a mixture of 
dark and bright states, in which the probability for the 
dark state is^ (I — e~ h °°)/2 and grows with temperature. 
As the lifetime of the dark state is infinite in this limit- 
ing case of parameter values the pure dephasing process 
effectively leads to a slowed down decay. 

The effect of thermal extension of the exciton lifetime 
is directly reflected in the dynamics of photoluminescence 
(PL) which may be observed in a time-resolved experi- 
ment. We calculate the PL signal by multiplying the 
occupations of the eigenstates |+) and |— ) by the corre- 
sponding decay rates T±. In order to account for ther- 
mally activated non-radiative processes, we follow Ref. 6 
and suppress the values obtained from our simulation by 
the factor exp(— T nr t), where r nr = T^t exp[— E a / {k^T)] 
is the rate of non-radiative exciton decay processes, at- 
tributed in Ref. [6j to exciton escape, which are here as- 
sumed to affect the two lowest exciton states in the same 
way. In this model, i? a is the activation energy for the 
non-radiative process and T^) is its overall strength. The 
PL decay time is then determined as the time at which 
the PL intensity (which decays not necessarily exponen- 
tially) is reduced by a factor 1/e. 

In Fig. 0Ja) we plot the decay time of the PL sig- 
nal from a single DQD structure after initial ultrafast 
preparation of the state (|f) + |2))/-\/2 (in fact, as we dis- 
cussed above, the initial condition is not very important 
if the phonon-induced thermalization is fast enough). We 
choose here a fixed value of A = 0.5 meV, E a — 95 meV 
(as estimated in Ref. @) and F^t = 8.16 ps _1 (half of the 
value in Ref. @, corresponding to a twice shorter radiative 
exciton life time assumed in our simulations). Varying 
the strength of the coupling results in a clear qualitative 
change of the temperature dependence of the PL decay 
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FIG. 4: (Color online) The decay time of the DQD photo- 
luminescence, (a) A single DQD with A = 0.5 meV and 
V = —0.2 meV (red solid line), —1 meV (black dotted 
line), —3 meV (blue dashed line), and —8 meV (green dash- 
dotted line), (b) An ensemble of DQDs with constant values 
of the coupling V as in (a), with the average energy mis- 
match A = 0, and with the standard deviation of A equal 
to a = 0.5 meV. (c) An ensemble of DQDs for A = and 
for three values the standard deviation of A: a = 1 meV 
(red solid line), 3 meV (blue dashed line), and 10 meV (green 
dash-dotted line). 



time. The shapes of the curves obtained from the simula- 
tions at lower temperatures reflect the rate of the radia- 
tive recombination, while at higher temperatures they are 
dominated by the non-radiative suppression. In general, 
the thermally increased quasi-equilibrium occupation of 
the higher, darker state leads to longer PL decay times, 
reduced at higher temperatures by the non-radiative con- 
tribution, which leads to the appearance of a maximum 
of the decay time. However, for a weak inter-dot cou- 
pling (red solid line in Fig. HJa)), the energy mismatch 
dominates and the two states are almost equally bright, 
that is, they have similar radiative decay rates. As a 
result, even at T = 0, when only the brighter state is 
occupied, the decrease of the life time is marginal and 
the resulting amplitude of the change in the PL decay 
time is small. This changes when the coupling is stronger 
(black dotted and blue dashed lines in Fig. Ufa)) so that 
the ground state has a considerably increased decay rate, 
which is manifested by a much shorter PL decay time at 
low temperatures. In a certain range of intermediate val- 
ues of V this results in a considerable growth of the PL 
decay time before the non-radiative effects become im- 
portant, which leads to the formation of a pronounced 
maximum. However, for higher coupling strengths the 
maximum disappears again. Now, the energy splitting 
between the two states, roughly equal to 2V, is so large 
that the change of equilibrium occupations is negligible 
in the temperature range in which the luminescence life 
time is governed by the radiative decay. 

As can be seen In Fig. EJb), a nearly identical de- 
pendence on temperature is obtained for an ensemble of 
DQDs in which the coupling V is again assumed fixed 
and the values of the energy mismatch A vary across the 
ensemble according to the Gaussian distribution with the 
mean value A = and a rather small standard deviation 
fj = 0.5 meV. The fact that the ensemble of DQDs is 
relatively homogeneous in terms of the energy mismatch 
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A is essential. As we show in Fig. @|c), if the standard 
deviation a becomes larger the shape of the tempera- 
ture dependence changes considerably. This is due to 
the fact, that for a 3> |V| the ensemble is dominated by 
DQDs with A 3> \V\ for which the collective effects are 
negligible. Hence, there is no decrease of the PL decay 
time at low temperatures for such ensembles. 

The temperature dependence of the PL decay time ob- 
tained from our model is qualitatively similar to that ob- 
served in the experiment 6 , where a reduction of the de- 
cay time by roughly a factor of two at low temperatures 
was observed and then a maximum was reached at about 
100 K, where the decay time equal to that observed in a 
reference sample with individual QDs was found. How- 
ever, as the above discussion shows, such a large ampli- 
tude of the variation of the PL decay time is only pos- 
sible for rather small values of V (in order to assure a 
small energy difference required for efficient thermaliza- 
tion). Therefore, our model is unable to quantitatively 
account for those experimental results obtained in a sys- 
tem containing much more strongly coupled DQDs, in 
which thcrmalization of occupations between the bright 
and dark states could not take place to a sufficient extent 
in the relevant temperature range. In view of the theory 
presented here, the origin of the strongly non-monotonic 
temperature dependence of the PL decay time in such 
a system remains an interesting issue. Based on our re- 
sults, we can only conclude that a thermally activated 
mechanism of breaking the coherent nature of the delo- 
calized exciton state, other than thermalization between 
the two states included in our model, is needed to account 
for this behavior. 



B. Transfer of coherence 

In this section we study the effect related to the res- 
onance between the biexciton shift Eb and the energy 
splitting between the single-exciton states. As discussed 
in Sec. IIII[ such a resonance affects the decay of inter- 
band polarizations in the DQD system, opening the path 
for a transfer of optical coherence from the biexciton 
transition to the single exciton transitions, according to 
Eq. (|9b]) . Here, we first discuss this effect in the absence 
of phonons. Then, we describe the phonon effects on this 
coherence transfer. Finally, we show that this dynamics 
of interband coherences is manifested by a line narrowing 
in nonlinear absorption spectra. 

In Fig. [51 we present the evolution of the interband co- 
herences (optical polarizations) for a single DQD in the 
situation when the single-exciton transition |— ) — > |0) is 
energetically close to the biexciton transition \B) — > \— ), 
without taking into account the carrier-phonon interac- 
tion. We plot the evolution of the excitonic polarization 
-Fx = \po-\, the biexcitonic polarization Px = \p~b\ and 
the total polarization in the spectral vicinity of the two 
transitions, defined as Ptotai = |po- + P-b\- Figs. [5Ja-c) 
show the dynamics for uncoupled dots (V = 0) at exact 
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FIG. 5: (Color online) The evolution of single-exciton, biexci- 
ton and total optical coherence for a DQD with 2£ = 3 meV 
in the absence of interaction with phonons. (a-c) Uncoupled 
dots with Eb = —3 meV, —3.01 meV, and —3.1 meV, respec- 
tively, (d-f) Coupled dots with 8 — — n/3 and with the same 
values of Eb as in (d-f), respectively. In all the panels the 
red solid line shows the total polarization, the green dashed 
line represents the biexciton polarization and the blue dotted 
line shows the single-exciton polarization. The evolution of 
the average exciton occupation (scaled down by a factor of 
3 to fit to the polarization amplitudes) is shown with gray 
dash-dotted lines for comparison. 



resonance, Eb — — 2 A = —3 meV and for two slightly 
off- resonant situations corresponding to \Eb\ — 2 A = 
0.01 meV and 0.1 meV. The initial state of the DQD is 
chosen to be (|-) + \B))/y/2. In all the three cases, the 
biexciton polarization (green dashed line) decays exactly 
in the same way (with the rate To + T_/2 = 3Fo/2 in 
this case). However, at resonance, its decay is accompa- 
nied by the appearance of a relatively long-living exciton 
polarization (blue dotted line). Hence, the decay time 
of the total coherence (red solid line) is considerably ex- 
tended. This effect is present only in a narrow energy 
range around the resonance. Already for a detuning of 
0.01 meV, the transfer of coherence to the single-exciton 
polarization is much smaller and the decay of the total 
polarization is only modulated by small oscillations. Fi- 
nally, at still larger detunings from the resonance, the 
coherence transfer vanishes completely and no increase 
of the polarization decay time is observed. As can be 
seen in Figs. [5j d-f), the evolution is very similar for cou- 
pled dots. Here all the time scales are slightly shorter 
since the involved single-exciton state is partly superra- 
diant. It should be noted that the strong variation of the 
polarization decay time at this particular resonance does 
not affect the evolution of exciton occupations, as repre- 
sented by the gray dash-dotted lines which show identical 
decay independent of the detuning. 

In Fig. [6l we plot the evolution of the same optical po- 
larizations as in Fig.[5Jbut taking into account the carrier- 
phonon coupling. In all the cases, the value of Eb is set at 
resonance (including compensation for phonon-induced 
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FIG. 6: (Color online) The evolution of single-exciton, biexci- 
ton and total optical coherence for a DQD with 2£ — 3 meV 
with the carrier-phonon coupling taken into account. The 
value of -Eb is taken at resonance, (a-c) Uncoupled dots at 
T = 0, 10 K and 20 K, respectively, (d-f) Coupled dots with 
6 = — 7r/3 at the same three temperatures. Line coding as in 
Fig. [5] 



energy shifts which are below 0.1 meV). Figs. a-c) cor- 
respond to uncoupled dots (V = 0) at the temperature 
T = K (a), 10 K (b), and 20 K (c). In this case, the 
evolution of polarizations is weakly affected by phonons 
due to the suppression of the phonon-induced real transi- 
tions between the two eigenstates. The only effect is the 
fast initial dephasing of the coherence^ which increases 
with temperature. In Figs. [BJd-f), we present the evo- 
lution of a DQD system with coupled dots (6 = — 7r/3, 
which corresponds to V « 1.3 meV and A w 0.75 meV). 
Here we see that the lifetime of the polarization is short- 
ened to some extent already at T = and becomes very 
short at T = 20 K. The amplitude of the single-exciton 
polarization induced by the polarization transfer also de- 
creases with temperature. The effect of pure dephasing 
is visible also here but it affects the evolution only over 
a very short period and leads to a small almost instan- 
taneous initial reduction of the polarization. The strong 
temperature dependence observed at longer times is due 
to real phonon-induced transitions between the single ex- 
citon eigenstates which is allowed in the case of coupled 
dots and is accompanied by a decay of coherences be- 
tween the single exciton states involved in these phonon- 
induced transitions and the other two states (vacuum 
and biexciton). Since our initial state involves only the 
lower-energy single exciton state the phonon scattering 
within the single exciton subspace induces the Marko- 
vian dephasing of the interband coherences with a rate 
which critically depends on the Bose occupation factor at 
the phonon frequency corresponding to the transition to 
the higher state, 71b(2£/7i). The major phonon-induced 
effect is therefore thermally activated. At higher temper- 
atures it reduces the lifetime of the polarizations from the 
usual time scales of radiative recombination to the values 
typical for phonon-related processes, that is, from almost 
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FIG. 7: (Color online) The pump-probe absorption spectra 
(see text) at T = 4 K (a), 10 K (b), and 20 K (c). Each plot 
shows three lines corresponding to exact resonance (red solid 
lines) and to a detuning of 0.05 meV and —0.05 meV from the 
resonance (blue dotted and green dashed lines, respectively). 
Here 2£ = 3 meV, 6 = -tt/3. 



a nanosecond to several picoseconds. 

While the analysis of the evolution of the individ- 
ual contributions to the polarization gives insight in 
the mechanisms underlying the system kinetics, these 
quantities are not directly accessed in an experiment. 
In order to study features related to the biexciton co- 
herences one usually resorts to nonlinear spectroscopy 
techniques2 ' 13 i 33 . For instance, the effects of dipole cou- 
plings between the dots have been thoroughly studied 
theoretically for the differential absorption measured in 
a pump-probe experiment^ 3 . Tracing the phonon effects 
on such spectra on the same level of exactness would re- 
quire including carrier-phonon coupling in the full model 
of Ref. 13 which would go far beyond the scope of this 
paper. Here, we only present a simplified model of a par- 
ticular pump-probe experiment in which the extended 
life times of the optical coherences are manifested in the 
absorption spectra in the most clear way. 

Let us consider a pump-probe experiment in which 
the spectrally resolved lowest-order response of a sin- 
gle pair of coupled dots to a probe pulse is studied af- 
ter a strong pump pulse (beyond the usual second or- 
der in the pump amplitude). The two pulses are spec- 
trally broad enough to excite all the relevant exciton 
and biexciton transitions in an identical way. The de- 
lay time between the pulses is long enough for the oc- 
cupations of the two single-exciton states to equilibrate 
but much shorter than the radiative life time. We as- 
sume that the area of the pump pulse a pU mp is such that 
sin 2 (a pump /2) = l/{2 - tanh[£ /(k B T)]}. This choice 
strongly suppresses the gain line corresponding to the 
emission from the ground exciton state leaving the biex- 
citon absorption line as the only feature in the relevant 
spectral range. We then calculate the probe absorption 
spectrum as the imaginary part of the Fourier transform 
of the total polarization obtained from the simulations 
with the initial state corresponding to an instantaneous 
excitation by the probe from the thermalized state pre- 
pared by the pump (in the linear order in the probe pulse 
area) . 
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The spectra obtained in this way, in the area around 
the transition energy from the lowest exciton state, are 
shown in Fig. [7] Here the energy is shown with re- 
spect to the average fundamental transition energy E — 
(ei +C2)/2. Each plot contains three spectra: one at res- 
onance (red solid line) and two detuned by ±0.05 meV. 
Each off-resonant spectrum contains the absorption line 
corresponding to the exciton-biexciton transition and a 
weak feature at the exciton-ground state transition en- 
ergy. The essential effect is the narrowing of the ab- 
sorption line at the resonance, which corresponds to the 
extended life time of the coherent polarization due to 
the coherence transfer to single-exciton polarization, as 
discussed above. This effect is quite strong at low tem- 
peratures but becomes less pronounced when the temper- 
ature increases and the line widths become dominated by 
phonon-related effects. 

V. CONCLUSION 

We have studied theoretically the kinetics of radiative 
recombination of excitons confined in systems consisting 
of two coupled quantum dots in the presence of carrier- 
phonon coupling. In general terms, this corresponds to 
spontaneous emission from non-identical coupled emit- 
ters undergoing an additional, possibly non-Markovian 
dephasing. We have shown that coupling to phonons 
strongly affects the dynamics of spontaneous emission by 
inducing transitions between the bright and dark states 
of the double-dot structure as well as by pure dephas- 
ing processes at short time scales. Both these phonon- 
related effects lead to strong temperature dependence of 
the carrier lifetime. We have shown that in the case of 



a binding ground state (negative coupling, typical e.g. 
for tunneling), the carrier lifetime increases with grow- 
ing temperature. For certain values of system parameters 
and in the presence of thermally activated non-radiative 
exciton decay processes, this can lead to a non-monotonic 
dependence of the exciton lifetime. 

We have also studied the process of resonant coherence 
transfer between the biexciton and exciton transitions 
and the effect of phonons on the kinetics of this process. 
We have shown that this effect can lead to considerable 
extension of the life time of coherent interband polar- 
ization in the double dot system, which is stable against 
phonon-induced decoherence in a relatively wide temper- 
ature range. We have presented a simplified model of a 
spectrally resolved pump-probe experiment which shows 
that the coherence transfer and the resulting extension of 
the polarization life time manifests itself by line narrow- 
ing in the nonlinear absorption spectrum of the system. 

Our results show that carrier-phonon interaction can 
qualitatively modify the optical response of double dot 
systems. Phonon-related effects turn out to play a much 
more important role in these systems than they do for 
individual dots. Our modeling can shed some light on 
the existing experimental results but also allowed us to 
predict new features that can be sought for in future 
experiments. 
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